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'. Abstract 

We analyze the hadronic and radiative decay modes of the recently observed 
_D s j(2317) meson, in the hypothesis that it can be identifield with the scalar 
sf = | + state of cs spectrum (D s0 ). The method is based on heavy quark 
symmetries and Vector Meson Dominance ansatz. We find that the hadronic 
^ . isospin violating mode D s q — > D s tt° is enhanced with respect to the radiative 

mode D s0 — > D*j. The estimated width of the meson is T(D s0 ) ~ 7 KeV. 
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1 Introduction 



The BaBar Collaboration has reported the observation of a narrow peak in the Dfir 
invariant mass distribution, corresponding to a state of mass 2.317 GeV [1]. The observed 
width is consistent with the resolution of the detector, thus T < 10 MeV. In the same 
analysis no significant signals are found in the D s '~f and -D s 77 mass distributions. The 
meson has been denoted as D sJ (23 17); the announcement has immediately prompted 
different interpretations [2, 3]. 

A possible quantum number assignment to _D s j(2317) is J p = + , as suggested by the 
angular distribution of the meson decay with respect to its direction in the e + — e~ center 
of mass frame. This assignment can identify the meson with the D s0 state in the spectrum 
of the cs system. Considering the masses of the other observed states belonging to the 
same system, Z) sl (2536) and _D s j (2573), the mass of the scalar D sQ meson was expected in 
the range 2.45 — 2.5 GeV, therefore ~ 150 MeV higher than the observed 2.317 GeV. A D s0 
meson with such a large mass would be above the threshold M^k = 2.359 GeV to strongly 
decay by S-wave Kaon emission to DK, with a consequent broad width. For a mass below 
the DK threshold the meson has to decay by different modes, namely the isospin-breaking 
D S 7T° mode observed by BaBar, or radiatively. The J p = + assignment excludes the final 
state -D s 7, due to angular momentum and parity conservation; indeed such a final state 
has not been observed. On the other hand, for a scalar cs meson the decay D s $ — > D*'-/ is 
allowed. However, no evidence is reported yet of the D s r y r j final state resulting from the 
decay chain D s0 — > D*j — > D s / -f-f. In order to confirm the identification of D s j(2317) with 
the scalar D s0 , one has at first to understand whether the decay modes of a scalar particle 
with mass of 2317 GeV can be predicted in agreement with the experimental findings 
presently available. In particular, the isospin violating decay to D s n° should proceed at a 
rate larger than the radiative mode D s0 — > .0*7, though not exceeding the experimental 
upper bound T < 10 MeV. This letter is devoted to such an issue. 

2 Mode D s0 D s n° 

In order to analyze the isospin violating transition D s0 D s tt° one can use a formal- 
ism that accounts for the heavy quark spin-flavour symmetries in hadrons containing a 
single heavy quark, and the chiral symmetry in the interaction with the octet of light 
pseudoscalar states. 

In the heavy quark limit, the heavy quark spin sq and the light degrees of freedom 
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total angular momentum sg are separately conserved. This allows to classify hadrons with 
a single heavy quark Q in terms of sg by collecting them in doublets the members of which 
only differ for the relative orientation of sq and sg. 

The doublets with J p = (0~, 1~) and J p = (0 + , 1 + ) (corresponding to sf = | and 



, respectively) can be described by the effective fields 

(1 + /) 
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where v is the four-velocity of the meson and a is a light quark flavour index. In particular 
in the charm sector the components of the field H a are = D^°, and D^f> (for a = 
1, 2, 3); analogously, the components of S a are P 0a = D®, Dq, D s0 and P[ a = Df, D'^~, D' sl . 

In terms of these fields it is possible to build up an effective Lagrange density describing 
the low energy interactions of heavy mesons with the pseudo Goldstone n, K and rj bosons 
[4, 5, 6, 7]: 
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In (3) H a and S a are defined as H a = ^H^ and S a = 7°S^7°; all the heavy field 
operators contain a factor y/M P and have dimension 3/2. The parameter A represents 
the mass splitting between positive and negative parity states. 

The 7r, K and 77 pseudo Goldstone bosons are included in the effective lagrangian (3) 

iM 

through the field £ = e f that represents a unitary matrix describing the pseudoscalar 
octet, with 
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and / ~ f n . In eq.(3) E = £ 2 , while the operators D and A are given by: 
D, ba = 5 ba d fl + V^ ba = 5 ba d fl + ^(^d^ + ^) ba 

Aba = Utfdpt-td^ 



(5) 
(6) 

The strong interactions between the heavy H a and S a mesons with the light pseu- 
doscalar mesons are thus governed, in the heavy quark limit, by three dimensionless 
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couplings: g, h and g' . In particular, h describes the coupling between a member of 
the H a doublet and one of the S a doublet to a light pseudoscalar meson, and is the one 
relevant for our discussion. 

Isospin violation enters in the low energy Lagrangian of n, K and r\ mesons through 
the mass term 

uf 2 

Cmass = ^j-Tr{£m 9 £ + (7) 
with m q the light quark mass matrix: 

/ rn u \ 
m q = m d . (8) 
V m s J 

In addition to the light meson mass terms, C mass contains an interaction term between 
7T° (7 = 1) and r] (I — 0) mesons: £ m i X ing — f nid ~^ lu it ti which vanishes in the limit 
m u = md- As in the case of D* — > D s 7r° studied in [8], the isospin mixing term can drive 
the D s0 — > -D S 7T° transition. 1 The amplitude A(D s0 — > D s tt°) is simply written in terms 
of A(D s o — > D s r]) obtained from (3), A(i] — > 7T°) from (7) and the rj propagator that puts 
the strange quark mass in the game. The resulting expression for the decay amplitude 
involves the coupling h and the suppression factor (m^ — m u )/(m s — md + m " ) accounting 
for isospin violation, so that the width T(D s0 — > D s n°) reads: 

rvn n 0\ 1 h 2 M Ds ( m d -m u x2 ^ 

r(^ — > D s n = — _ (1 + [ ^ F )kol . (9) 

As for h, the result of QCD sum rule analyses of various heavy-light quark current 
correlators is \h\ = 0.6 ± 0.2 [6]. Using the central value, together with the factor 
(m d - m u )/(m s - HMip) ~ _i_ [9] anc i f = f w = 132 MeV we obtain: 

T(D s0 ^D s n°)~QKeV. (10) 

Eq.(9) can receive SU(3) F corrections: a hint on their size comes from the use of / = 
= 171 MeV instead of f n in (9), which gives T{D s0 -> D s tt°) ~ 4 AW. On the other 
hand, we neglect corrections related to the finite charm quark mass. 

The analogous calculation for D* — > D s n° involves the coupling g in (3). Since h and 
g have similar sizes (0.3 < g < 0.6), it turns out that the transitions D s0 — > D s tc° is en- 
hanced with respect to D* — > -D s 7r° essentially due to kinematics, being ^j^^^oyp 
3 x 10 2 . 
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1 Electromagnetic contributions to £) s0 — > -DsTr are expected to be suppressed with respect to the 
strong interaction mechanism considered here. 



3 Radiative D s q 



D*7 decay 



Let us now turn to D s0 — > -D*7, the amplitude of which has the form: 

A(D s0 -> L>: 7 ) = erf [(e* • V *)(p ■ k) - ( V * ■ p)(e* ■ k)] , (11) 

where p is the D s0 momentum, e the D* polarization vector, and k and 77 the photon 
momentum and polarization. The corresponding decay rate is: 

F(D s0 -> D* sl ) = a\d\ 2 \k\ 3 . (12) 

The parameter d gets contributions from the photon couplings to the light quark part 
e s s7^s and to the heavy quark part e c c7 M c of the electromagnetic current, e s and e c being 
strange and charm quark charges in units of e. Its general structure is: 

d = dW + dW = |L + ^-, (13) 

where A a (a = c, s) have dimension of a mass. Such a structure is already known from 
the constituent quark model. In the case of Ml heavy meson transitions, an analogous 
structure predicts a relative suppression of the radiative rate of the charged D* mesons 
with respect to the neutral one [10, 11, 12, 13], suppression that has been experimentally 
confirmed [14]. From (12,13) one could expect a small width for the transition D s0 — > -D*7, 
to be compared to the hadronic width D s0 — > D s ir° which is suppressed as well. 

In order to determine the amplitude of D s q — > D*^ we follow a method based again 
on the use of heavy quark symmetries, together with the vector meson dominance (VMD) 
ansatz [11, 13]. We first consider the coupling of the photon to the heavy quark part of 
the e.m. current. The matrix element (D*(v', e)\c^ffj,c\D s o(v)) (v, v' meson four-velocities) 
can be computed in the heavy quark limit, matching the QCD cry^c current onto the 
corresponding HQET expression [15]: 

jHQet = ~ hv[v ^ + *(# _ *Q ) + -^—a^{T + T) + . . ]h v (14) 
ZnriQ Zttiq 

where h v is the effective field of the heavy quark. For transitions involving D s0 and D*, 
and for v — v' (v • v' — 1), the matrix element of J^ QET vanishes. The consequence is 
that d^ is proportional to the inverse heavy quark mass vtlq and presents a suppression 
factor since in the physical case v • v' — (m 2 DsQ + m 2 D *) /2m Ds0 m D * = 1.004. Therefore, we 
neglect d^ in (13). 
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To evaluate the coupling of the photon to the light quark part of the electromagnetic 
current we invoke the VMD ansatz and consider the contribution of the intermediate 
0(1020): 

(15) 

with k 2 = and (O|s7 M s|0(/c, ei)) = M^f^e^. The experimental value of is = 
23AMeV. The matrix element (D*(v' : e)<j>(h, €i)\D s0 (v )) describes the strong interaction 
of a light vector meson (0) with two heavy mesons (D* and D s0 ). It can also be obtained 
through a low energy lagrangian in which the heavy fields H a and S a are coupled, this 
time, to the octet of light vector mesons. 2 The Lagrange density is set up using the hidden 
gauge symmetry method [5], with the light vector mesons collected in a 3 x 3 matrix p M 
analogous to M in (4). The lagrangian 3 reads as [16]: 

C! = tpTr{S a H b a x »V Xu (p) ba } + h.c. , (16) 

with V Xv {p) = d x pu - d u p x + \p\,p v ] and p x = i^Px, 9v being fixed to g v = 5.8 by the 
KSRF rule [17]. The coupling fx in (16) is constrained to p = —0.13 ± 0.05 GeV^ 1 by the 
analysis of the D — * K* semileptonic transitions induced by the axial weak current [16]. 
The resulting expression for j- is: 

1 _ A ,,9v [M^ U (17) 



The parameters are obtained from independent channels; we use their central values. 
The numerical result for the radiative width: 

T(D s0 -> D* sl ) ~ 1 KeV (18) 

shows that the hadronic D s0 — > D*tt° transition is more probable than the radiative mode 
D s0 — > -D*7- In particular, if we assume that the two modes essentially saturate the D s0 
width, we have 

T(D s0 )~7 KeV (19) 

and 

B(D s0 -> D s 7r°) ~ 0.85 (20) 
B(D s0 ^D* sl ) ~ 0.15 



2 The standard oj s — uj mixing is assumed, resulting in a pure ss structure for <j>. 
3 The role of other possible structures in the effective lagrangian contributing to radiative decays is 
discussed in [13]. 
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at odds with the case of the D* meson, where the radiative mode dominates the decay 
rate. 

The same conclusion concerning the hierarchy of D s0 — > D s ir° versus D s0 — > D*^ is 
reached in [3] using the quark model. Since our calculation is based on a different method, 
the sf = \ and sf = | + doublets being treated as uncorrected multiplets, we find the 
agreement noticeable. 



4 Conclusions and perspectives 

We have found that the observed narrow width and the enhancement of the D s ir° decay 
mode are compatible with the identification of D sJ (2317) with the scalar state belonging 
to the sf = | + doublet of the cs spectrum. However, this conclusion does not avoid 
other questions raised by the BaBar observation, one being the low mass of the state. We 
believe that such a particular issue requires additional investigations. A second point is 
that the radiative mode, although suppressed, is not negligible, and should be observed 
at a level typically represented by the ratios in (20). 

The quantum number assignment has two main and rather straightforward conse- 
quences. The first one is the existence of the axial vector partner D' sl belonging to the 
same spin doublet sf = | + . Even in the case where the hyperfine splittings between 
positive and negative parity states are similar: M D > si — M Ds0 ~ M D * — M Ds , this meson 
is below the D*K threshold. Therefore, its hadronic decay to D*ir°, at the rate 

F(D> - D**°) - ^ Md ° ( md ~ mu \ 2 \2 + (Mk + M ^ " M '°) 2 1 



x (l + ^)b";o| 3 ^r( J D s0 ^ J D s 7r ) (21) 

would produce a narrow peak in the D*n° mass distribution. The confirmation of such a 
state, the existence of which is suggested by the analysis of the Ds^yn mass distribution 
[1], will support the interpretation. 

The second consequence concerns the doublet of scalar and axial vector mesons in 
the bs spectrum. Since the mass splitting between B and D states is similar to the 
corresponding mass splitting between B s and D s states, such mesons should be below 
the BK and B*K thresholds, thus producing narrow peaks in B s ir° and B*n° mass 
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distributions, with rates resulting from expressions analogous to (9)-(21). 



Note added. When this work was completed, the CLEO Collaboration announced the 
observation of a narrow resonance with mass 2.46 GeV in the D* + 7i° final state and the 
confirmation of D s j(2317) [18]. Moreover, a theoretical analysis based on the quark model 
was posted on the Los Alamos arXive, with the same conclusions presented here [19]. 

Acnowledgments. We acknowledge partial support from the EC Contract No. HPRN- 
CT-2002-00311 (EURIDICE). 
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